Long-lived oscillons from asymmetric bubbles: existence and stability 
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The possibility that extremely long-lived, time-dependent, and localized field configurations ("os- 
cillons") arise during the collapse of asymmetrical bubbles in 2+1 dimensional (j)'^ models is investi- 
gated. It is found that oscillons can develop from a large spectrum of elliptically deformed bubbles. 
Moreover, we provide numerical evidence that such oscillons are: a) circularly symmetric; and b) 
linearly stable against small arbitrary radial and angular perturbations. The latter is based on a 
dynamical approach designed to investigate the stability of nonintegrable time-dependent configu- 
rations that is capable of probing slowly-growing instabilities not seen through the usual "spectral" 
method. 
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I. INTRODUCTION 



The existence of classical field configurations exhibit- 
ing soliton-like properties is not only an interesting con- 
sequence of nonlinear effects in field theory, but also an 
important ingredient in the understanding of nonpertur- 
bative effects in particle physics Q. Apart from one- 
dimensional kinks, examples of these structures in 3-f 1 
dimensions include nontopological solitons (NTS) and 
Q balls 1^: in both cases, choosing a simple harmonic 
time-dependence for the scalar field allows one to ob- 
tain a static solution of the field equations, describing a 
spherically-symmetric configuration which, for a range of 
parameters, may be the lowest energy state. Such con- 
figurations may be found in extensions of the standard 
model, supersymmetric or not, as has been suggested 
recently Q. They may be sufficiently stable as to al- 
low for a quantization procedure and form a legitimate 
bound state (see e.g. Q). NTSs and Q-balls have also 
been of great interest to applications of particle physics 
to the early universe, often being proposed as possible 
candidates for dark matter (5|. A time-dependent, long- 
lived, and localized configuration in 3-fl dimensional 
4>'^ scalar field theory was re-discovered and thoroughly 
studied by one of us and collaborators some years ago 
1^, 0]. It was shown that these configurations, named 
in Ref. |^ "oscillons," naturally arise from collapsing 
unstable spherically-symmetric bubbles in models with 
symmetric and asymmetric double- well potentials, being 
mainly characterized by a rapid oscillation of the field 
at the bubble's core. Their relevant feature is that, al- 



beit not strictly periodic, they possess a very long life- 
time, of order 10^ — IO^to"^ (?i = c = 1 in this work). 
Oscillons may be thought of as the higher-dimensional 
cousins of one-dimensional breather states found from 
kink-antikink bound states Just as kink-antikink 

pairs may be thermally or quantum-mechanically nucle- 
ated through nonperturbative processes, so may oscillon 
states, although here the calculation must be done in real 
and not Euclidean time. More recently, Gleiser and Sorn- 
borger investigated whether oscillons are present in 2+1 
dimensions, finding not only that they do exist, but also 



that their lifetime is at least of order 10 



Moti- 



* Electronic address: [artur. adib@brown.edu ; Address after Septem- 
ber 2002: Box 1843, Brown University. Providence RI 02912, USA. 
t Electronic address: g'leiserQlogos. dartmouth.edu 
t Electronic address: harIos@fisica. ufc.br 



vated by this result, in the present paper we investigate 
two important related questions: first, if oscillons still 
appear during the collapse of asymmetric - as opposed 
to symmetric - initial configurations; and, second, if they 
are stable against small angular and radial fluctuations. 
Since this implies that we will be dealing with a higher- 
dimensional parameter space, we restrict ourselves here 
to 2-1-1 dimensional oscillons. Apart from being of in- 
terest in their own right, we expect that our results will 
be indicative of the behavior of 3+1 dimensional oscil- 
lons. We also note that it should be quite easy to build 
oscillons from more complicated field theories, including 
interactions between the "oscillon" field and other scalar 
or fermionic fields. The robustness of these configura- 
tions, as demonstrated here, should provide enough mo- 
tivation for a careful search of such generalized oscillons 
(and possibly more realistic) in the near future. 

The paper is organized as follows. By means of a nu- 
merical scheme suited to tackle long-lived configurations 
(described in the Appendix), in Sec. || we show that 
oscillons quickly appear during the collapse of most el- 
liptically deformed bubbles and, moreover, that they are 
all circularly symmetric and extremely long-lived, leav- 
ing no trace of the initial asymmetry. This suggests 
that oscillons can be understood as attractors in field- 
configuration space, ordered spatio-temporal structures 
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that emerge during the nonhnear evolution of a wide va- 
riety of initial configurations. In fact, in 2+1 dimensions, 
the attractor basin is quite deep, as was initially hinted in 
Ref. joj and will be further shown here. We then move on 
to study, through a dynamical approach, whether these 
symmetric configurations are stable against small asym- 
metric perturbations, finding no indication of spectral in- 
stability (Sec. [II). We conclude in Sec. IV summarizing 
our results and pointing out future avenues of research. 



II. OSCILLONS FROM ASYMMETRICAL 
BUBBLES 

The Lagrangian density for our 2-1-1 dimensional scalar 
field theory is: 



(1) 



with fi — 0,1,2. We introduce dimensionless variables 
by rescaling the coordinates and the field as — X ^TTL 

and (j)' — (f)^/X/m (henceforth we drop the primes). The 
energy and the equation of motion are 



El 



and 



-{dcj^/dtf- 



I? 



(2) 



(3) 



respectively. So far, all previous studies have obtained os- 
cillons from symmetric initial configurations, with either 
thick or thin walls (Gaussian or tanh profiles, respec- 
tively). We will restrict our investigation to Gaussian 
initial profiles, as these proved to be the most interest- 
ing in 2-1-1 dimensions (cf. ||^). For convenience, we will 
follow Ref. j|] and restrict the initial field configuration 
to interpolate between the two minima of the potential. 
Of course, one could select different values for the initial 
value of the field at the core [(j){t = 0,r — 0,6)]: as it 
was argued in Ref. , as long as the value of the field at 
the core probes the nonlinearity of the potential, and the 
initial configuration has an energy above the "plateau" 
energy (the energy of the oscillon configuration) , oscillons 
are bound to appear. 

The asymmetry in the initial field configuration is in- 
troduced by means of an elliptical deformation: 



where: 



0(p,6',O) = -2exp 



Rl 



1, 



1 — cos^ ' 



(4) 



(5) 



with e : [0, 1) the bubble eccentricity, Rq the bubble "ra- 
dius" and p, 6 polar coordinates (notice that this expres- 
sion reduces to the usual symmetric ansatz when e = 0). 
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FIG. 1: Logarithm of the "effective radial dispersion" 
ARcf{(t = 0) vs. eccentricity e for the ansatz (^. Ai?cff 
is clearly a good measure of asymmetry, increasing always 
monotonically with e. 
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FIG. 2: Oscillon energy in a shell of radius Rs — 5Ro for 
-Ro = 2.0 and 0.1 < £ < 0.8, from bottom to top (only integer 
multiples of e = 0.1 are shown). 



We note in passing that a similar parameterization was 
adopted in the study of eccentric pulsons in the sine- 
Gordon theory ||l^ . 

In order to measure the asymmetry of the field con- 
figuration, we take advantage of the effective radius i?off 
defined in Ref. and introduce an "effective radial dis- 
persion" A_Reff (see below) . The effective radius could be 
recast in the form: 



i?eff(i) = {R{0,t))g, 



(6) 



where {■)g = (27r)-i d^ is an angular average and 



3 




FIG. 3: Oscillon energy in a shell of radius Rs — 5Ro, now 
for Ro = 3.0 and 0.1 < e < 0.8, with steps Ae = 0.05. Dashed 
lines indicate unstable configurations. 



FIG. 4: Oscillon energy for Ro — 4.0 and different e, with 
steps Ae = 0.05 (dashed lines are for unstable configurations). 



R{6, t) is a ^-dependent effective radius defined by (com- 
pare Eq. (28) of Ref. 0): 



R{9,t) 



Cdp [l(V0)^ + y(0)] 
(2^)-i^'"d0/J^=dpp[i(V^)2 + y 



(7) 



with Rg the "shell radius" (i.e. the radius within which 
we integrate all the quantities relative to the osciUon). 
This notation aUows us to construct angular averages and 
dispersions in a manner analogous to time averages. We 
therefore define the effective radial dispersion as the rela- 
tive root-mean-square deviation from the average radius: 



Ai?eff(i) 



^(RHe,t))e-{R{e,t)yi 

{Ri9,t))g 



(8) 



As shown in Figure the above quantity is indeed a 
good measure of asymmetry, i.e. it increases monoton- 
ically with e. ft is approximately independent of the 
bubble size, being the limiting case Ai?cft = an indica- 
tion of a symmetrical state (though not exactly zero on 
a lattice due to its finite resolution). We now turn to the 
presentation of the main numerical experiments obtained 
by solving Eq. (^ with the initially asymmetrical bub- 
bles (^ for eccentricities e ranging from 0.1 to 0.8 and Rq 
from 2.0 to 5.0. This investigation can easily be extended 
to greater values of Rq, although this would require much 
longer computational times without generating results of 
much physical interest. [The computational time Ctimo 
is proportional to L^, which in turn is proportional to 
(i?o)^, see Appendix]. 

Figures |^|[ show the time evolution of the total en- 
ergy within a shell of radius Rg = 5i?o surrounding the 
initial configuration for different values of Rq and e. It is 
seen that, in general, initially asymmetric bubbles tend to 
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FIG. 5: Shell energy for bubbles of radius Ro = 5.0 and steps 
Ae = 0.05 (dashed lines for unstable cases). 



decay into coherent field configurations with an approx- 
imately constant energy plateau similar to those found 
in Ref. which focused on the evolution of symmet- 
ric configurations. With the help of the effective radial 
dispersion Ai?off, we can investigate whether these con- 
figurations correspond to "excited" states of an oscillon 
(i.e., non-spherically symmetric configurations analogous 
to an excited state of a hydrogen atom for £ ^ 0) or if 
the bubble asymmetry is completely lost and the system 
decays into a "ground" (i.e. symmetric) state. 

In Fig. II we show the time evolution of Ai?off for 
the Rq ~ 4.0 case. It is clearly seen that the initially 
asymmetric bubble decays into a AR^s = 0, symmetric 
configuration, after a brief asymmetric pulsation. A sim- 
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FIG. 6: A7?eff for Ro = 4.0 and 0.1 < e < 0.6, from bottom 
to top (only integer multiples of e = 0.1 are shown). Notice 
how the initial asymmetric bubble decays into a circularly 
symmetric configuration. 
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FIG. 7: Collapse of ARes in Figure | using AR^gie, t) = 
e°'f{t) for a ~ 2.05. Note that the time domain is restricted 
to t < 250, since Ai?eff — > for t > 250 (causing the log scale 
to be ill behaved). 



ilar pattern was observed for all initial radii investigated 
here, suggesting that whenever an oscillon stage is set, 
the resulting configuration is circularly symmetric. The 
peaks in this figure also suggest that Ai?cff might follow a 
scaling law for different values of e. Indeed, a collapse of 
these curves using the scaling Ai?off(e, t) — e"f{t), where 
a is a real constant and f{t) is a function of time only, 
is shown in Figure |^. 

As the reader must have noticed, an intriguing feature 
of these results is the presence of some "instability win- 
dows" for some values of e. These can be observed here 
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FIG. 8: Phase portrait between t — and t — 1000 at p = 
of a bubble entering an oscillon stage for Rq = 2.0 and £ = 
(top) and e = 0.4 (bottom): the symmetric configuration to 
which asymmetric bubbles decay is an oscillon. 



in the cases Rq — 3.0 and Rq — 5.0. Thus, oscillons 
do not always appear as the asymmetric configurations 
decays. A finer investigation of the parameter space for 
the elliptical deformations, generalizing what was done in 
great detail for spherically-symmetric 3-1-1 dimensional 
oscillons ||ll|, will quite possibly reveal a very rich and 
detailed substructure of stable and unstable windows. It 
is important to stress that once AR^s (cf. fig. ||), 
the field does settle into an oscillon, as the phase space 
portrait of Fig. ^ exemplifies. This justifies our earlier 
claim that oscillons are attractors in field-configuration 
space. 

A crucial step not yet studied is the stability analy- 
sis against small but arbitrary asymmetric perturbations. 
The oscillon stability with regard to these perturbations 
is fundamental for the computation of quantum correc- 
tions around the classical solution ^ and is the subject 
of the next session. 



III. STABILITY ANALYSIS 

Our task now is to investigate whether symmetric oscil- 
lons are stable against small radial and angular perturba- 
tions S(t){p,0,t), i.e. to probe the "linear stability" of the 
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oscillon. We are unaware of any previous study where the 
stabihty of a field configuration with a time-dependent 
ampHtude has been tested against small fluctuations. 
So far, stability investigations have been restricted to 
either time- independent configurations, such as bounce 
solutions []l5| , or to configurations with a linear time- 
dependence in the phase, such as Q-balls [^). The sta- 
bility analysis in these two situations is greatly simplified 
by the fact that the dynamical equation dictating the be- 
havior of the fluctuations is separable into its spatial and 
temporal parts; the resulting problem reduces to finding 
the eigenvalues of a time-independent operator (an 
alternative approach based on the so-called "Bogomolnyi 
bound" jl^ can also be applied in the time-independent 
case, see e.g. [|l3|). The existence of at least one negative 
eigenvalue (and thus of a complex eigenfrequency a;„) 
signals the presence of an exponentially-growing instabil- 
ity [Q (a well-known example is the so-called "bounce" 
solution, which we will investigate further below [^ ). 
The present problem, however, is not amenable to such 
treatment due to the anharmonic time-dependence of the 
oscillon; we must consider both the space- and time- 
dependence of the background field, making the stability 
analysis of oscillon- type configurations considerably more 
challenging both analytically and numerically, as we will 
now discuss. 

In order to appreciate these difficulties, let us write (in 
polar coordinates) the linearized equation of motion that 
follows from Eq. (||) through the substitution (/)(p, 9, t) — > 
(j)o{p,t) ~\- S(f>{p,6,t), where (po is the symmetric oscillon 
solution and S4> <C 4>o is the perturbation, i.e.. 



5^ 



dp- 



■5(1)- 



1 d , 

pop 



p2 96*2 



<50-(3(/.2-l)<50. (9) 



Here one might be tempted to separate the variables as 
Sep = R{p)Q{6)T{t). However, the resulting equations 
show that one cannot get rid of the simultaneous radial 
and temporal dependence of the background configura- 
tion, (j}o{p,t) This situation should be contrasted 
to the usual case where (po is a time-independent solu- 
tion, e.g. the bounce, or to the case where the time- 
dependence of 00 is in a phase factor exp{iivt), and thus 
immediately eliminated in the full equation of motion . 
A considerable simplification can nevertheless be accom- 
plished by writing d(j) = ^{p, t)Q{9), isolating at least the 
angular part of the problem. Performing such substitu- 
tion gives the pair of equations. 



„ 1 , 

P 



Uo{p,t) 



$ 0, 



and 



(10) 



(11) 



with overdots and primes indicating time and radial 
derivatives, respectively. Here Uo{p,t) = 30o — 1 and £ 
is a separation constant. The solutions for O are triv- 
ial, viz. Q cx exp(±i^0), and by requiring S(j) to be 



single- valued we have £ = 0,1, 2, .... Our original (2-1-1)- 
dimensional problem, Eq. (^, reduces therefore to solv- 
ing the above (l-l-l)-dimensional one, Eq. (p^, since the 
time-dependence is present only in $. 

Our goal is to probe the linear stability of the oscil- 
lon by solving Eq. (10) for arbitrary initial conditions. 
The strategy is to monitor the time evolution of the per- 
turbations d(j), which should grow without bounds in the 
case of a linearly unstable configuration |jl^, |l^ . An ob- 
vious limitation of this approach is that it is impossible 
to scan all initial values of perturbations, viz. 6(j){p,0) 
with £ = 0,1,2,.... The method thus is only indica- 
tive of stability, not being able to provide conclusive 
proof. The more thorough the search, the more one is 
guaranteed to show stability, at least against most types 
of perturbation. This unavoidable limitation should be 
contrasted with the simpler case for time-independent 
background configurations based on a harmonic decom- 
position T{t) cx exp{iujnt) (see e.g. (|, |4|), where the 
existence of exponentially unstable modes is clearly re- 
lated to imaginary eigenvalues. However, we would like 
to point out a limitation of the spectral method that is 
often overlooked. By restricting the analysis to an expo- 
nential time-dependence, as in T{t) above, one can ob- 
tain only spectral instabilities of a configuration, leaving 
aside other possible forms of instabilities, for example, 
linear (or power-law) ones. In other words, a system 
that is spectrally stable may still be unstable against 
slower growing modes |8j. Since we are here es- 
sentially watching the full time-dependence of dcj), we 
should be able to detect any sort of instability by ob- 
serving its long-time behavior, although in practice the 
infinite-time limit or a complete scan of possible fluctu- 
ations cannot be achieved numerically. Fortunately, we 
shall shortly see that typical spectral instabilities (such 
as that of the bounce) do not require a long-time inte- 
gration or a very wide search, being therefore bound to 
be observed through our method. Before we do so, it is 
worth testing the reliability of the numerical implemen- 
tation itself. 



A. Linear test 

The first step is to compare the numerical solution 
of Eq. (|l^) with a closed-form analytical one in or- 
der to prepare and test our numerical implementation, 
since the singular behavior at the origin requires a care- 
ful treatment. This can be done most easily by set- 
ting Uo{p,t) = 0, in which case Eq. ( pO| ) becomes 
linear and separable, with $ = R{p)T{t). This gives 
T cx exp(±jwi), where w is a separation constant, and 
the equation for R, 



R" + -R' 



u;^-^l]R = 0, 



(12) 



which we recognize as Bessel's equation. By requiring 
regularity at the origin and $(/0, 0) = 0, the solution can 



6 




The bounce solution 4>b{p) is the 0(2)-synimetric static 
configuration that solves the equation 
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FIG. 9: Semi-log plot of the maximum amplitude ^max ver- 
sus time for the bounce solution (showing here three different 
initial conditions). The slope denotes the common unstable 
mode eigenvalue, uin- 



be written as 

^t{p-,t)^ I dLuf{uj)Je{Lup)cos{ujt), (13) 



1 d(f> 
P dp 



with the asymmetric potential 



(16) 



(17) 



In order to detect the instability, we solved Eq. (|l^) with 
Uo{p) = f3—24>i,{p)+i<pl{p) and various initial conditions, 
sweeping the lattice at every time step to find the maxi- 
mum value of the perturbation, ^max- In Fig- ^ we show 
our results for /? = 0.011 [the initial conditions are Eq. 
(|l|) with n = 0,l and m = 2, and Eq. (|^) for m = 2, 
both with £ = 0] ; one can clearly identify the exponential 
growth of S(p even at early times t < 100. Also shown is 
the slope of the curve, which should match the unstable 
eigenvalue w„ obtained with the usual spectral method. 
[We have attempted to obtain such eigenvalue by solv- 
ing numerically the associated Schrodinger-like equation. 
However, in two spatial dimensions the severe singularity 
at the origin causes a numerical instability which we were 
unable to control even with sophisticated methods [^ . 
Since this is not the focus of this paper, we will leave this 
question aside.] 



where /(w) is determined by the initial condition. Choos- 
ing $(p, 0) = AJp[ap) we have 



/(w) ^Auo I dppJt{ap)Ji{up) = A6{lu ~ a), (14) 
Jo 

and therefore 



$f(p, i) = A dw S{lj ~ a)Jf{ujp) cos{ujt) 



AJi(ap) cos(at). 



(15) 



The above solution maintains its shape but oscillates har- 
monically with period 27r/a. We have verified that our 
numerical implementation reproduces correctly this ana- 
lytical solution for various values of a and £. 



B. The bounce 

As a first application of our method we investigate the 
stability of the so-called "bounce" solution [|l^ , which is 
guaranteed to be spectrally unstable in any dimension 
greater than (1+1) due to Derrick's theorem [^. (In 
fact, Coleman has showed that only one negative eigen- 
value exists |2^). Should our method be reliable, the so- 
lution S4> for the case where 0o is a bounce solution will 
grow exponentially at late times, indicating the presence 
of an unstable mode. 



C. The oscillon 

We are now ready to apply our method to the stability 
of the oscillon, which was obtained here by solving Eq. 
(^) with the symmetrical version (e = 0) of the ansatz 
Eq. (^ . We have essentially followed the same procedure 
described above for the bounce, but now evolving both 
Uq and $ in Eq. (0). Since the dimensionality of the 
configuration space is infinite, we chose arbitrarily the 
initial profiles of the fluctuations ^{p,to), with the only 
constraint that they should vanish at p ^ oo to ensure 
localization around the oscillon. The time was chosen 
to be about 200, since that is roughly when the initial 
bubbles have just decayed into an oscillon (cf. Figs. 
Some examples of the initial configurations investigated 
here are 



<^>i{p,to) = a;"exp 



11 

Hq 

, , , . / 2mrp\ 
'^i{p,to) = sm ( I cxp 

$Kp,to) = JKp)exp ( - — 



(18) 
(19) 
(20) 



for various integer values of n, m and £ (namely, n — 
0,1,2, TO = 2, 3 and € = 0,1, 2, 3). 

In Figure ^ we show a typical outcome of our search. 
In all cases investigated we found that the fluctuations 
^max are bounded from above, as one would expect from 
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1000 



FIG. 10: A typical outcome of the linear stability analysis 
of the oscillon, showing here $max vs. time for ^ = 0, n = 
and m = 2 in Eq. (flq). The radius of the initial bubble is 
Ro = 2. 



a linearly stable configuration. We conclude that if, in- 
deed, there are any unstable modes, they are sufficiently 
slow-growing to justify the use of the oscillon as a stable 
bound state. [We have integrated the linearized equa- 
tions of motion up to i = 10"*; see also discussion below]. 
Note that the large amplitude of $, e.g. ^ 1000 in Fig. 
|l0| , does not mean that the condition S<j) ^ 0o is violated: 
since the resulting equation is linear, any solution can al- 
ways be rescaled without changing its shape by choosing 
a different constant prefactor for the initial conditions. 

Although not as systematic and transparent as the in- 
vestigation above, another approach to check the sta- 
bility of oscillons is to superimpose the perturbation to 
the full (2-|-l)-dimensional oscillon dynamics discussed in 
Sec. ||. One can then probe the oscillon stability sim- 
ply by checking the persistence of the energy plateau: 
if the added energy from the perturbations is radiated 
away, the oscillon is stable. Due to the dimensionality of 
the problem, the numerical treatment is quite more chal- 
lenging than the one use above within the linear method. 
Nevertheless, we have investigated the stability of oscil- 
lons against superimposed fluctuations for similar initial 
conditions. In Fig. ^ we present the outcome of a par- 
ticular choice of initial condition for three different initial 
radii. The results are consistent with the previous sta- 
bility analysis, as can be seen by the persistence of the 
energy plateau. 

On the basis of our extensive search with many dif- 
ferent initial conditions and long integration times, we 
find it very unlikely that an exponentially-growing mode 
exists. If it does, it would be either very small and/or 
related to a very "rare" excitation; the oscillon config- 
uration would still be stable for large times and could 
be considered a legitimate (or at least a very long-lived) 



FIG. 11: An example of the full dynamics of the oscillon 
when subject to a perturbation of the form in Eq. ( |l8| ) for 
n — and m = 2 a.t t = 300. The maximum amplitude 
of the initial perturbation - here the prefactor in Eq. ( p^ - 
was constrained to either 0.01 or 0.05, although the results do 
not change appreciably even for 0.1. The stability is evident 
through the persistence of the energy plateau. 



bound state in semi-classical quantization. 



IV. CONCLUSIONS 

We have investigated, in 2-1-1 dimensions, two key 
questions concerning the properties of oscillons - time- 
dependent, localized field configurations that emerge dur- 
ing the deterministic evolution of 0^ models. First, 
we have shown that initially asymmetric configurations 
evolve, for a wide range of elliptic deformations, into sym- 
metric oscillons states. Thus, oscillons are not just par- 
ticular to symmetric initial states. This result led us 
to propose that oscillons are attractors in field configu- 
ration space, with a very deep attractor basin, at least 
in 24-1 dimensions. Second, we have shown that oscil- 
lons are stable against a wide range of asymmetric small 
perturbations. This result was obtained by two distinct 
approaches, one solving the linearized equation for the 
perturbations and the other by superimposing the per- 
turbations on the oscillons and evolving the perturbed 
configurations with the full equation of motion. Clearly, 
both methods are restricted to the choice of initial fiuc- 
tuations. However, after an extensive search, we were 
unable to find any unstable fiuctuation with either ap- 
proach. To the best of our knowledge, this is the first 
dynamical investigation of the stability of explicitly time- 
dependent scalar field configurations. We expect that 
both of these results will carry on to 3-1-1 dimensions, 
although probably the attractor basin will be shallower 
in this case. 



These results suggest the importance time-dependent 
spatio-temporal structures may have in a wide range of 
physical systems, from condensed matter to early uni- 
verse cosmology. Although we have restricted our study 
to simple (j)"^ models, we expect, as suggested in Refs. 
[0, H], that oscillons will be present whenever there is a 
bifurcation instability related to the negative curvature 
of the nonlinear potential. Oscillons will emerge in a 
wide variety of dynamical systems, possibly representing 
a bottleneck to equipartition of energy, thus delaying the 
approach to equilibrium. 

One possible arena for oscillons in early universe cos- 
mology is during the reheating supposed to occur after 
inflation. Oscillons may be thermally nucleated with a 
probability proportional to exp[— E'osc/T'], where E^sc is 
the energy of the oscillon configuration. They will act as 
"entropy sinks" , confining several degrees of freedom to 
an ordered state, delaying the thermalization of the uni- 
verse. Eventually, when they decay into radiation, they 
will dump more entropy to the early universe, possibly 
changing the final reheating temperature. 

Finally, it would be interesting to compute the spec- 
trum of quantum fiuctuations around oscillon states, 
to investigate their effect on oscillon stability. Time- 
dependent bound states may have much to add to 
our knowledge of nonperturbative quantum field theory, 
which has traditionally focused on time- independent con- 
figurations, such as instantons. 
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Put briefly, the leapfrog equations with the damping 
method read: 



1 + r7ijAt/2 



(A.l) 



where superscripts (subscripts) denote temporal (spatial) 
indices, rji^j is the damping function of Ref. P] and 
V'{(l)) is a first partial derivative of the potential with 
respect to the field. The second spatial derivatives in 
the Laplacian are discretized with a fourth-order scheme 

(to wit, dxx4>i,] ~ [16(</'j+lj +4>i-l^j) - 4>i+2,j ~ 4>i-2,j - 

30(^ij]/12Aa;^ and analogously for dyy(j)ij), which gives 
an energy conservation of one part in 10^ for Ax = 0.1 
and At — 0.06 (and, of course, with rjij = 0). A bet- 
ter energy conservation could be obtained with smaller 
Ax or At, but this comes with a high price tag since, 
as remarked below, the computational time is inversely 
proportional to both Ax^ and At. Despite this fact, with 
the above parameters we were able to reproduce quite ac- 
curately the results of Gleiser and Sornborger ||]. Even 
though the damping method is already a major improve- 
ment over more naive methods (such as huge lattices or 
even moving boundary conditions), for the problem at 
hand it is still demanding. As an example, for small os- 
cillons of radius Rq « 2.0, the required lattice of radial 
dimension R « 200 adopted in Ref. ^ (and thus L « 400 
in our square grid, where L is the lattice edge), would al- 
ready demand a total of ~ 10^ sites for Aa; = 0.1, as 
opposed to the N ~ 10"^ used in the latter reference. An- 
other aggravating fact comes from the large integration 
times T involved in such problems [notice that the re- 
quired computational time for this problem goes roughly 
asCtimo ~ {T/At)N = {t / At){L / Ax)^]. We note in pass- 
ing that there has been some effort to find a more natural 
and efficient discretization for the (f/^ theory which might 
reduce significantly the computational time of such prob- 
lems Motivated by this possibility, two of us have 
recently investigated these lattices and have found that, 
unfortunately, they are of limited practical use even for 
simple dynamical problems ||2^ . It was seen, however, 
that if the above scheme is supplied with the boundary 
conditions of Ref. [E4l : 



APPENDIX: NUMERICAL METHOD 

The integration scheme adopted here is a standard 
leapfrog algorithm, which ensures second-order preci- 
sion in time p2| (the spatial discretization used is 
fourth-order). We have adopted the "adiabatic dam; 
ing method" (or simply the damping method) of Ref. 
together with Higdon's first-order boundary conditions 
^4| . Their combined use turned out to be very effec- 
tive and of easy implementation, allowing us to tackle 
this otherwise demanding numerical problem with cur- 
rent workstations. 



(A.2) 



where a is either x or y, then a significantly smaller lat- 
tice could be used, resulting in an energy error smaller 
than (or equal to) the error due to numerical energy fluc- 
tuations. [These first-order "absorbing boundary condi- 
tions" were obtained for the rather simple (linear) wave 
equation. We expect, however, that the damping intro- 
duced before the boundaries could reduce the amplitude 
of the outgoing waves such that Eq. (^) is effectively 
linearized in t hat region, and thus that the boundary 
condition (A.2) becomes applicable]. With regard to the 
example in the previous paragraph, we have found that 
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the required lattice with this "mixed" method needs only 
L 140 (in contrast to the former L k, 400), such that 
N (and thus Ctimc) is roughly one order of magnitude 
smaller than the previous one (this trend is also found for 
greater i?o)- For the sake of completeness, we quote here 
the parameters of the damping method used through- 
out our simulations (we use the same functional form 
for 7y(p) as Ref. k = 0.005 (damping constant). 



Po = 10i?o (initial radius of the damping) and pi = 50 
(damping length), these latter two being defined such 
that L = 2(po + pi)- 

We expect that the method adopted here might be 
useful not only in higher dimensional systems (the two 
methods above do not really make any dimensional re- 
quirement), but also in other finite-domain problems not 
necessarily related to oscillons. 
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